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^ ABSTRACT 

We present VLT/VIMOS IFU observations of an occulting galaxy pair previously 

,— ^ discovered in HST observations. The foreground galaxy is a low-inclination spiral 

^^ disk, which causes clear attenuation features seen agains the bright bulge and disk of 

r \ the background galaxy. We find redshifts of z = 0.064 ± 0.003 and z=0.065 for the 

^^ foreground and background galaxy respectively. This relatively small difference does 

,^ not rule out gravitational interaction between the two galaxies. Emission line ratios 

Mh point to a star-forming, not AGN-dominated foreground galaxy. 

^ We fit the Cardelli, Clayton & Mathis (CCM) extinction law to the spectra of 

\^ individual fibres to derive slope (i?y) and normalization {Ay). The normalization 

"^ agrees with the HST attenuation map and the slope is lower than the Milky Way 

Ci relation {Rv < 3.1), which is likely linked to the spatial sampling of the disk. We 

'~~' speculate that the values of Ry point to either coherent ISM structures in the disk 

, larger than usual (^ 9 kpc) or higher starting values of Ry, indicative of recent 

►^ processing of the dust. 
ly-v The foreground galaxy is a low stellar mass spiral {M^, ~ 3 x lO^M©) with a 

y—i high dust content (Mdust ^ 0.5 x 10^ Mq). The dust disk geometry visible in the 

\^ HST image would explain the observed SED properties of smaller galaxies: a lower 

C^ mean dust temperature, a high dust-to-stcllar mass ratio but relatively little optical 

_il attenuation. Ongoing efforts to find occulting pairs with a small foreground galaxies 

(^— ^ will show how common this geometry is. 

, Key words: techniques: imaging spectroscopy techniques: spectroscopic galaxies: 

. . distances and redshifts galaxies: dwarf galaxies: individual: 2MASXJ00482185-2507365 

t> galaxies: ISM galaxies: spiral 
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1 INTRODUCTION perature gradient in the dusty ISM (Bendo at al. 2010; Smith 

et al. 2010; Engelbracht et al. 2010; Pohlen et al. 2010; Foyle 

The distribution and properties of interstellar dust in spi- , , „mr,\ ^ i_i • ii • i /n j ^ i nnir,N 

, , r , et al. 2012), notably m the spiral arms (Bendo et al. 2010) 

ral galaxies is still the dominant unknown m models ot the i i i /t^ h , . . i o^if>\ r , , ■ i 

r , , and bulge (Engelbracht et al. 2010) ol nearby galaxies, and 

spectral energy distribution (SED) or spiral galaxies (e.g., ,, ^ i-j_i totv-t i i ^ /r-. i ^ i 

"■^ o V o ' the power spectrum or the ISM, gas and dust (Combes et al. 

Popescu et al. 2000; Popescu & Tufls 2005; Pierini et al. r,^ir,\ tt ^^ n j_ i ^-n i i i i 

„. „ , . '„ „ ^. ,. , , ^. ,. 2012 . However, the smallest scales can still only be probed 

2004; Baes et al. 2003, 2010; Bianchi et al. 2000; Bianchi '^ ' ^ f 



2007, 2008; Driver et al. 2008; Jonsson et al. 2010; Holw- 



in some of the most nearby galaxies, e.g., Andromeda, the 

Magellanic Clouds or the Milky Way itself. In the case of 
erda et al. 2012a). Recently, the focus has shifted to how ,, i • rr ^ o rn i , tt i i n 

/ / n-m-N smaller galaxies, opiizer 6pace lelescope and Herscnel space 

the dust IS distributed, and its impact on SED models and ^, , ,, . . . i , r ■ i- • i i i 

Observatory results are restricted to a lew individual galax- 

distance measurements (e.g., SNIa, Holwerda 2008). ^^^ ^^ ^^ ^.^^ ^^ ^^ 2006, 2012; Hermelo et al. 2013), and 

Knowledge of the dust geometry in external spira ^^^^^ ^^^^ ^^^^^^^ ^^^^ p^p^^^^ ^^ ^^ 2002; Grossi et al. 

galaxies has seen great improvements with the Herschel r,r>io -n ^ i r,oir,\ mi i ^- ■ ^i j_- ^ 

r 1 , 2010; Bourne et al. 2012). thus, observations m the optical 

Space Observatory. This new tacihty has mapped the tem- i^ • i j_ i ^i • j::ii i ^- ^ 

" J j-i- Qj. ultra-violet wavelength regimes can hll an observational 

gap to estimate the dust geometry in spiral disks or dwarf 

galaxies independently from sub-mm observations 
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To map dust extinction or attenuation^, one needs a 
known background source of light. In the case of a spiral 
galaxy overlapping a more distant galaxy, the more distant 
galaxy provides such a background light source. Estimating 
dust attenuation and mass from differential photometry in 
occulting pairs of galaxies was first proposed by White & 
Keel (1992). Their technique was then applied to all known 
pairs using ground-based optical images (Andredakis & van 
der Kruit 1992; Berlind et al. 1997; Domingue et al. 1999; 
White et al. 2000), spectroscopy (Domingue et al. 2000), 
and later space-based HST images (Keel & White 2001a,b; 
Elmegreen et al. 2001; Holwerda et al. 2009). These results 
agree with distant galaxy counts through spiral disks in HST 
images, an alternative approach (Holwerda 2005; Holwerda 
et al. 2005a,b,c,e,d, 2012b). More recently, additional pairs 
were found in the SDSS spectroscopic catalog (86 pairs in 
Holwerda et al. 2007), and in an effort in the Galaxy Zoo 
project (Lintott et al. 2008); the current count is 1993 pairs 
(Keel et al. 2013). Holwerda et al. (2009) (hereafter H09) 
have recently found a superb example of an overlapping pair 
in Hubble Space Telescope (HST) data (Figure 1). 

This system is one of the best examples of an occulting 
pair: it combines an optimal viewing configuration (face- 
on foreground disk and partial overlap) and a bright back- 
ground galaxy. This pair was serendipitously observed with 
the Advanced Camera for Surveys (ACS) in observations of 
the outskirts and halo of NGC 253 (Bailin et al. 2011) as 
part of the ACS Nearby Galaxies Survey Treasury (ANGST, 
Dalcanton et al. 2009). Our analysis in H09 showed an ex- 
tended dusty disk in the foreground galaxy, extending well 
beyond its optical size. The average reddening-attenuation 
relation appears to be close to the Milky Way Extinction 
Law {Rv = {Ab — Av)/Av = 3.1), although with significant 
scatter (Figure 2 in H09). This measurement was based on 
the broad-band FST filters {F475W, F606W and F8I4W) 
and averaged over the entire overlap region. The attenuation 
is largely clustered in a series of spiral structures (Figure 2 
in H09). Such an extended dusty disk poses several ques- 
tions for current models of ISM in spirals: how did so much 
dust end up (or form) so far outside the optical disk? Is this 
common or is this particular galax;y an anomaly? Is this pair 
truly a chance occultation or interacting in some way? 

To address these questions, we present here 
VLT/VIMOS-IFU observations of this pair with the 
specific aim to obtain the redshift difference between both 
galaxies, as well as a measurement of the attenuation 
curves in a distant spiral disk. This paper is organized 
as follows: section 2 describes the VIMOS observations, 
sections 3 and 4 discuss the observed redshifts and emission 
line class of the foreground galaxy respectively, section 
5 discusses two observational effects in determining the 
observed attenuation curve, section 6 present our analysis 
of the VIMOS fibres, and section 7 is our discussion of these 
results. 



^ Because the actual contribution of scattered liglit is unknown, 
we will use the term dust attenuation throughout the paper, ex- 
cept where discussing the extinction-reddening relations, i.e., the 
observed extinction curves and law. 



2 VLT/VIMOS OBSERVATIONS 

The Visible MultiObject Spectrograph (VIMOS Le Fevre 
et al. 2003) on the VLT/UT3 can operate in three different 
modes; either as an imager, a Multi-Object Spectrograph, 
or as an Integral Field Unit (IFU), which is the mode used 
for this study. The IFU consists of 4 quadrants of 1600 fibres 
each, feeding four different 2kx4k CDDs. Each quadrant is 
made by four sets (pseudo-slits) of 400 fibres. We asked for 
O'.'67/fibre (27" x 27" field-of-view) with arcsecond seeing or 
better, to resolve as much detail of the overlap region as 
possible while covering the entire object. The medium spec- 
tral range mode (MR) is 4900-10150 A with R=720 spectral 
resolution with the G4T5 filter. 

The galaxy pair is listed in the 2MASS (Klein- 
mann et al. 1994) catalog under 2MASXJ00482185-2507365 
(RA=00:48:21.8, DEC=-25:07:37) with a major axis diam- 
eter of 0^5. During ESO period 84A (on 2009-12-18), 1.2 
hours of service mode observations with VLT/VIMOS in 
IFU/MR mode were obtained of the target pair (program 
384.B-0059A). Seeing measurements were reported between 
~ 0'.'7— 0'.'8, and the airmass was less than 1.02. The point- 
ing successfully covered most of the pair (Figure 1). 

Our observing strategy was based on a very tentative 
identification of the foreground galaxy's redshift as z=0.03 
based on an HIPASS (Meyer et al. 2004) radio spectrum. 
Thus, we opted for medium spectral resolution to maximize 
throughput and obtain the best signal-to-noise on the even- 
tual attenuation curve. The spectral resolution of the MR 
spectral mode is not high enough to obtain kinematic infor- 
mation for either galaxy. 

The VIMOS pipeline delivers high-level science prod- 
ucts for observations in service mode. We use the flux- and 
wavelength-calibrated images, and the extracted fibre spec- 
tra^. The standard pipeline performs bias subtraction and 
extracts fibre-spectra, after solving for wavelength and dis- 
persion using sky-lines, and resampling them to a constant 
wavelength step per pixel (2.6 A per pixel). Individual spec- 
tra are corrected for fibre transmission and re-ordered in a 
spectral image for each quadrant with each fibre-spectrum 
in a row. The corresponding sky positions can be found in 
the position table ifutablesHR.tar (from the ESO website), 
which we combined into the field-of-view image in Figure 1. 
We assigned WCS sky coordinates to this fibre map based 
on the assumption that the position of the brightest fibres 
correspond to the nuclei of both galaxies. Small-scale cor- 
rections such as pixel-to-pixel response ("flat-fielding") and 
fringing correction are not executed in the VIMOS pipeline 
(Izzo et al. 2004) or subsequent improvements (Scodeggio 
et al. 2005). In addition to the pipeline reduction we flat- 
fielded and subtracted the sky emission as discussed below. 



2.1 Flat Fielding 

As we are interested in the extinction curve by taking the 
ratio between fibres, the fibre spectra would not need be flat- 
fielded, provided the fibre response corrections are identical 
for both fibres. However, this is certainly not the case in 



■^ see also http://lf s.wikidot .com/vlmos (Westmoquette et al. 
2009) and the ESO's website http://www.eso.org/. 
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HST/ACS (F606W) 



VIMOS R=720 




Figure 1. The HST/F608W iraage (left) and the VLT/V/MOS reconstructed image (right) of the occulting galaxy pair (top two panels) 
and the sky-subtracted spectra of both nuclei (bottom panel). The blue and red spectrum correspond to the fibres on the centers of 
the foreground {FG, blue circle) and background galaxy {BG, red circle), respectively. The dotted spectrum is the sky contribution 
estimated from several "empty" fibres in the fourth quadrant of the VIOS IFU. This sky contribution was subtracted from both FG 
and BG spectra. We analyze two apertures of the overlap region; one in which the foreground galaxy contributes a significant amount 
of flux (I, green box) and one for which we assume no foreground disk emission (II, black box). The flux contributions to Aperture I are 
estimated from the corresponding flbre positions for the foreground galaxy {F' , blue box) and background galaxy (C, red box). Optical 
depths are calculated from equation 1. Optical depth estimates for Aperture II (black box) are estimated from the ratio of flux in the 
aperture and the corresponding fibres on the other side of the background galaxy (orange box) and equation 2. 



VIMOS. Small-scale corrections add both noise and compli- 
cate the fit to the slope of the extinction curve. Flatfielding is 
done on an individual fibre basis using the quadrant-specific 
flat-field files (MXFF). The lamp spectrum for a fibre is re- 
trieved, and a third order polynomial is fit and divided out 
to remove the lamp's black body spectrum. This normal- 
ized residual is the response function for fiat-fielding and 
the spectrum is divided by this function. 

Zanichelli et al. (2005) present a VIMOS graphical 
pipeline and point out that to correct for overall fibre re- 
sponse, one can use the emission in the 5577 A sky line, 
which should be uniform across the field-of-view. We checked 
using this line and the relative total fibre responses are uni- 
form within a few percent. 



Table 1. The identified lines in the spectrum of the foreground 
galaxy's nucleus (blue circle and spectrum in Figure 1). 



Line 


(A) 


(A) 


FWHM 

(A) 


EW 


z 


Line 

Strength 

(counts) 


[Hell 


5876 


6380.6 


11.75 


2.10 


0.0859 


0.15 


[OI] 
fNIIl 


6300 
6548 


6670.5 
7182.5 


5.09 
4.71 


1.64 
4.22 


0.0588 
0.0671 


0.10 
0.22 


Ha 


6563 


6987.4 


4.75 


5.73 


0.0647 


0.33 


fNIIl 


6583 


6968.1 


4.44 


-0.09 


0.0642 


-0.01 


[SII] 


6731 


7008.4 


6.33 


4.62 


0.0646 


0.24 



2.2 Sky Estimate and Subtraction 

After flat-fielding of each fibre as explained in section 2.1, the 
nigh sky spectrum needs to be subtracted. Zanichelli et al. 
(2005) point out that the VIMOS instrument does not have 
dedicated fibres to estimate the sky background. However, 
in cases where the objects do not cover the entire field-of 
view, such as ours, they advocate the use of the median of 
those fibres with lowest total intensity as a good measure 
for the sky background. This is what we have done for each 
quadrant; we computed the average of the group of fibres in 



each quadrant, which are in the lower range of the summed 
fibre flux values. 



3 INTERACTING OR NOT? THE REDSHIFTS 
OF BOTH PAIR MEMBERS 

We estimate the redshifts of both galaxies based on the fibre 
spectra of their respective nuclei (the blue and red circles 
in Figure 1, respectively). Figure 1 also shows the usable 
spectrum from these fibres (4900-7300 A). 
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We identified five emission lines in the foreground 
galaxy's spectrum (although the [He] redshift identification 
is uncertain, Figure 1). We fit a Gaussian profile to each 
line and used the central wavelength to estimate a redshift 
for the foreground galaxy (Table 1). The mean redshift is 
z = 0.0639 ± 0.0027 for the foreground galaxy (excluding 
the [He] line). 

Two absorption feature redshifts for the combined ob- 
ject (2MASXJ00482185-2507365) are reported in the lit- 
erature; z=0.063641±0.0005 (Ratcliffe et al. 1998) and 
z=0.064011±0. 000150 (6dF survey, Jones et al. 2009). These 
are essentially redshift estimates for the background galaxy 
because it dominates the flux. Both previous studies identify 
the spectrum as that of a bulge or elliptical, dominated by 
old stellar spectra. 

In the spectrum of the background galaxy's nucleus (red 
circle and aperture in Figure 1), we identified two absorption 
features (Mg and Na, indicated in Figure 1). Their position, 
compared to template spectra for a Bulge or SO (Calzetti 
et al. 1994; Kinney et al. 1996), confirm the previous red- 
shift estimates {z = 0.06495±0.0007). In the outskirts of the 
background galaxy's disk, some faint star formation is ap- 
parent in the HST image. However, it is not strong enough 
for a certain emission line identification in the relevant fibre- 
spectra. 

The difference in redshift between the pair members is 
very small; Az = 0.001 (Av^yg = 300 km/s), or only some 4 
Mpc, assuming Hubble Flow (Ho — 73 km/s/Mpc). 

This separation is still substantial but the relative veloc- 
ities do not completely preclude an interaction in this pair's 
history or future. Indeed, both pair members could easily 
be at the same distance and hence strongly interacting. In 
our opinion, this is unlikely, because both galaxies appear 
otherwise undisturbed and the dusty structures appear to 
form a symmetric disk around the foreground spiral, not a 
tidal feature. 



4 STARBURST OR AGN? LINE STRENGTHS 
IN THE FOREGROUND DWARF GALAXY 

If we compare the line strength in Table 1, we find ratios 
of -0.18 and -0.14 for log([NII]/Ha) and log([SII]/Ha), re- 
spectively. Following the Baldwin et al. (1981); Kewley et al. 
(2006) classification (BPT diagram), these line ratios iden- 
tify the foreground galaxy as most likely an HII (star form- 
ing) object, not an AGN or LINER. If the foreground galaxy 
has an AGN, it is a fairly weak one (we lack H/3 informa- 
tion to make a complete classification). If both galaxies are 
indeed in close proximity, the current star-formation in the 
foreground galaxy may well be tidally triggered. 



5 OBSERVATIONAL EFFECTS ON THE 
OBSERVED EXTINCTION CURVE 

There are two observational effects to consider when com- 
paring the observed relation between attenuation and wave- 
length in the overlap region: scattered light into the line-of- 
sight and the effects of physical sampling of the foreground 
ISM. Asymmetry in either galaxy introduces uncertainty in 



the extinction measures. Instrumental effects are discussed 
in section 6.3. 



5.1 Scattering 

Because the galaxies are not a point sources but extended 
objects, the relation between attenuation and wavelength 
may include effects of light scattered into the line-of-sight by 
dust clouds in the foreground disk. This effect is treated in 
detail in the appendix of White et al. (2000) . There are two 
sources of scattered light, the foreground and background 
galaxy. In Aperture I, the foreground galaxy contribution is 
removed (see below), and in Aperture II, this contribution 
is much lower than the background galaxy's. This leaves 
the background galaxy as the main source of scattered light 
in our measurement. Scatter would affect our measurement 
most where the background source is brightest, i.e., close 
to the background galaxy's nucleus, specifically Aperture II 
(see below). 

White et al. (2000) examine the contribution of scat- 
tered light into a line-of-sight as a function of angular and 
radial separation of the two galaxies. If the galaxies are suffi- 
ciently far apart, the scattering angle becomes small enough 
not to matter. They express the scattering percentage as a 
function of the ratio between radial extent and separation of 
both disks. The foreground galaxy disk has a radius of ~12 
kpc and the separation implied by the respective redshifts is 
4 Mpc. The ratio is thus several hundred which results in a 
scattering effect well under a percent of the measured flux. 



5.2 Sampling 

In H09, we already noted the effect of spatial sampling on the 
observed extinction law. Assuming the background galaxy's 
redshift as an upper limit, we explored the effects of resam- 
pling the image on the attenuation measure. Now that the 
actual distance of the foreground disk is known, we can ex- 
plore, using the HST attenuation maps as a start, predict 
the spatially averaged extinction curve that will be observed 
in the VIMOS spectra. 

Figure 2 shows the reddening {E{B — V)) versus atten- 
uation (Av) in the HST attenuation maps (Figures 8 and 
9) for the native resolution and the resampled case where 
the HST maps are spatially smoothed to the 0'.'8 seeing and 
resampled at 0'.'67 pixels of the VIMOS observations. The 
spatially smoothed and resampled maps show a very gray 
observed extinction curve {Rv « 3.1), averaged over the 
whole HST overlap region at Ry = 1.3. This implies that 
the notion of measuring the extinction curve as a function 
of position in the foreground disk of the foreground galaxy 
may be too ambitious for the actual sampling and resolution 
of the observations. 

We adopt the term extinction curves and law for the 
remainder of the paper with the caveat that these measure- 
ments are influenced by sampling of the disk, with some 
scattered light from the background galaxy as a secondary 
effect in the very brightest sections of the background galaxy. 
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Figure 2. The reddening, E{B — V) = Ag — Ay versus attenuation in F606W, Ay determined for each pixel in the attenuation map 
in Figure 9. For the optical depth maps at the native HST resolution a relation of Ry = 2.86 is retrieved (fit to the positive values of 
Ay). To compare with the VIMOS observations, we spatially smooth and resample the HST images to O'.'S resolution with 0'.'67 pixels 
(right panel) to explore how the distance and seeing affect the extinction law that would be observed with 'VIMOS resolution. Spatially 
smoothed and resampled images yield lower values of Ry, i.e., a grayer extinction law. The right panel shows the extinction law one 
would expect in the VIMOS observations based on the resampled ilf5T images (FWHM=0'.'8, pixelscale=0'.'67); Ry = 1.29. 



5.3 Galaxy Asymmetry 

The central assumption in our analysis of the HST and VI- 
MOS observations (below) is that both galaxies are rota- 
tionally symmetric. This assumption of symmetry must es- 
pecially hold for the background galaxy. For this reason, we 
would prefer an elliptical galaxy as a background. However, 
the background galaxy in this case appears a very regular 
spiral without much structure from star-formation which of- 
ten breaks symmetry on small scales. In principle, the dark 
dust structure visible in the HST image and the attenuation 
evident in the VIMOS data (Figures 1, 7 and 9) could be part 
of the background galaxy but it would uncharacteristically 
break symmetry for the background galaxy. And secondly, 
the general morphology of the dusty structures strongly sug- 
gests that they belong to the foreground disk (orientation, 
implied winding angle) and not the background galaxy. 



6 OBSERVED EXTINCTION CURVES 

We analyze two parts of the overlap region: one where we 
have an estimate of the foreground galaxy's contribution 
(Aperture I in Figure 1), and one where we assume there is 
little flux from the foreground galaxy (Aperture II in Figure 
1). Counterpart fibres for both foreground and background 
galaxies are retrieved from the blue and red or orange Aper- 
tures in Figure 1 respectively. We compute the attenuation 
as a function of wavelength from: 



A{\) 



-1.086 X In 



[i?(A)e-(^)+f(A)]-[F'(A)] 
[B'(A)] 



(1) 



where B and B' are the contributions from the background 
galaxy (overlap and counterpart respectively) and F and F' 
the contributions from the foreground galaxy (overlap and 
counterpart). The terms in square brackets correspond to 
an observable spectrum: overlap [Be^^ + F], foreground [F'] 



and background [B'] counterparts. The implicit assumptions 
are that both galaxies are intrinsically rotationally symmet- 
ric and that the foreground counterpart fibre (F') has little 
fiux contributed by the background galaxy due to the the ex- 
ponential drop in background galaxy surface brightness with 
radius. As explained above, equation 1 ignores the effect 
of scatter and the effect of mixing transparent and opaque 
lines-of-sight. 



6.1 Aperture I: Foreground Stellar Disk 

In Aperture I (green box in Figure 1), the flux contribution 
from the foreground galaxy (F) is not negligible and we use 
equation 1 to estimate the attenuation in the foreground 
disk ([Be~^ -I- F], foreground [F'] and background [B'] are 
the green, blue and red box in Figure 1 respectively). The re- 
sulting extinction curves from the raw spectra are uncertain 
and noisy (Figure 3), as the flux contribution from the fore- 
ground galaxy is flrst subtracted and the background galaxy 
is relatively faint, compounding the issue of individual pixel 
signal-to-noise. A viable extinction curve cannot be derived 
from the individual, raw spectra. 

There are two ways to mitigate this, either by spectrally 
smoothing the individual flbre spectra, sacriflcing spectral 
resolution, or by combining multiple fibres, sacrificing spa- 
tial information. First we choose the first approach; Figure 3 
shows the result of smoothing one of the fibres spectrally by 
21 pixels (300 A), to improve the signal-to-noise in the ex- 
tinction curve. However, only a very noisy extinction curve 
can be obtained. Even if we omit subtracting the foreground 
galaxy spectrum, no realistic extinction curve can be re- 
trieved. 

Thus, we can combine the spectra from the two least 
noisy fibres in Aperture I into a single average extinction 
curve (Figure 5, green line), but even this "mean" extinction 
curve holds very little information. 
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Figure 3. An example of a raw (left) and spectrally smoothed (right) fibre spectrum in Aperture I. The spectra were Hanning smoothed 
by 21 pixels (54. 6A) to improve signal to noise. After aggressive spectral smoothing, a noisy extinction curve can be recovered. However, 
the CCM fit to this curve is problematic for both the slope and normalization. 
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Figure 4. An example of a raw (left) and spectrally smoothed (right) fibre spectrum in Aperture II. The spectrally smoothed spectrum 
was convolved with a 21 pixel Hanning profile. The CCM fits to both extinction curves agree well {Ry = 0.9, Ay = 0.9). The normalization 
is close to the values derived from the HST attenuation map in the FdO&W filter (Figure 9). The fit to the slope {Rv) is lower than 
expected from the HST broad-band filters. 



6.2 Aperture II: Foreground Outer Disk 

In Aperture II (black box in Figure 1), our analysis is 
simpler, because we can ignore the flux contribution from 
the foreground galaxy, as the background galaxy dominates 
[F « B). Aperture II encompasses the overlap region we 
analyzed in H09 (Figure 8). We estimate the optical depth 
as a function of wavelength in these fibres from: 



A{X) = -1.086 X In ( i|^^ 



(2) 



Individual extinction curves for the 15 fibres of Aper- 
ture II are shown in Figure B2. We note that in the majority 
of fibres in this Aperture, there is enough s/n even in the 
unsmoothed spectra to generate a (noisy) extinction curve. 
The one exception where a CCM fit was problematic (Table 
Al and Figure B2), fibre F15, is the one closest to Aper- 



ture I (Figure Al), where it proved impractical to infer an 
extinction curve. We can attempt to improve the s/n in the 
extinction curve in the same two ways as above: by spectrally 
smoothing individual fibre spectra in the spectral dimension 
and/or by combining spectra, in effect averaging spatially. 

Figure Bl illustrates the effect of spectral smoothing an 
individual fibre spectrum. While the noise is reduced, the fit 
to the extinction curve (see below) remains the same (but 
with lower x^). The mean extinction curve from the aver- 
aged spectra in all of Aperture II (and the corresponding 
background galaxy aperture) is shown in Figure 5. By sacri- 
ficing the spatial information, the noise can also be brought 
down. One thing is immediately evident: the spatially av- 
eraged extinction curve is very flat, i.e., a gray extinction 
curve, in line with what one would expect for an average of 
the whole aperture from §5.2. 
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Aperture I (F3+4) 

Aperture II 

Aperture I CCM (J?v.= 1.3, Ay= 2.1, x' =416.9) 

Aperture II CCM («,= 1.9, Ay= 0.5, x' = 22) 




6000 6500 

Wavelength (A) 



Figure 5. The two mean extinction curves for botli apertures; 
green line (the two best fibres in the green box in Figure 1) and 
black line (Aperture II in Figure 1). We show fits of the CCM re- 
lation to these two extinction curves. Aperture I shows extremely 
high Ay values and an unexpectedly high value of Ry . The mean 
values for Aperture II are in line with expectations from the HST 
aperture, both in normalization (Ay) and slope (Ry). 




1.5 2.0 2.5 

Normalization A^ 



Figure 6. The Ay and Ry values with the spread due to different 
sky-subtraction. Each quadrant's sky estimate was subtracted in 
turn from the overlap and comparison spectrum and the rms of 
the resulting CCM fits is shown as the uncertainty. 



6.3 Uncertainties and Biases 

In section 2, we point out that we supplement the VI- 
MOS pipeline with sky-subtraction and a fibre small-scale 
sensitivity correction. In the analysis above, we compare 
the fluxes different fibres to each other. Differences in the 
fibre-to-fibre transmission, sky-subtraction or wavelength re- 
sponse (flat-fielding) could affect the resulting fibre extinc- 
tion curves, both the normalization (Ay) and the slope 

The spectral response for individual fibres can affect 
both the slope (Ry) and normalization (Ay) of the derived 
extinction curve. The relative correction for overall fibre re- 
sponse, done in the pipeline, for both the fibre in the over- 
lap region and the comparison one(s) in either galaxy affects 
the normalization of the extinction curve. In addition, the 
small scale, pixel-to-pixel correction affect the overall fit. 
For this reason, we supplemented the pipeline with individ- 
ual, normalized, fibre flat-fields. Flat-fielding predominantly 
changes the noise and spectral slope of the fibre spectra. If 
both the fibre in the overlap (either aperture) and the back- 
ground galaxy fibre were affected in the same manner, it 
would only make the CCM fit somewhat noisier but there 
are subtle differences in the slopes of the VIMOS fibre flat- 
fields, which could cause a bias in the slope, Ry as well. 

Sky-subtraction has the potential to change both the 
normalization and slope of the extinction curve but only 
if there are spatial variations in the sky background. The 
field-of-view of the VIMOS/IFU is sufliciently small that 
one would not expect any real variation in the sky emission, 
a fact we used to check the individual fibre transmission. 
However, each quadrant of the IFU has a slightly different 
throughput in the optics and for this reason we compute sky 
estimates for each quadrant. 

Fringing is most evident in the red part of these VIMOS 
spectra (the VIMOS received an upgrade in 2010 to coun- 
teract this feature as well as flexure). In principle, this does 



not affect our extinction curve estimate except that it makes 
the spectra and resulting extinction curve more noisy. 

Differences in the sky-subtraction or transmission be- 
tween fibres would affect the resulting extinction curve. This 
relative response calibration is well taken care of in the VI- 
MOS pipeline with the relative transmission file. Zanichelli 
et al. (2005) point out that the VIMOS instrument does 
not have dedicated fibres to estimate the sky background. 
However, in cases where the objects do not cover the en- 
tire field-of view they advocate the use of the median some 
fibres with low total intensity as a good measure for the 
sky background. This is what we have done for all our fibre 
spectra; subtracted an appropriate sky background for each 
IFU quadrant (e.g., Q4 in Figure 1). To estimate the effect 
of sky subtraction, we repeated the CCM fit with the sky 
estimate for each quadrant subtracted from the overlap fibre 
and the corresponding background fibre in Aperture II. The 
rms spread in both values, the normalization and slope is 
shown in Figure 6. Typical rms values in the normalization 
are Ay < 1 but with larger uncertainties in the slope, Ry. 



6.4 Mapping the Extinction Curve 

We fit the Cardelli, Clayton and Mathis (CCM) extinction 
curve to each of the fibre extinction curves (really attenua- 
tion curves, see section 5) in Aperture II. Figure 4 illustrates 
the effects of spectrally smoothing the fibre-spectra; the fit 
improves but generally the results did not change for each 
fibre (see also Figures A2 and Bl). Only for the extreme 
value of Ry is there a marked difference between the spec- 
trally smoothed and the unsmoothed fit results. We there- 
fore map the fit results to the unsmoothed fibre extinction 
curves in Figure 7 (also listed in Table Al). 

The slope of the extinction curve (Rv) is between 1.5 
and 4.2 in most realistic cases (Figure 7, right), with the 
majority of the values lower than the canonical Ry = 3.1 
but often higher than we found by averaging the whole HST 
aperture (Figure 8). A plausible explanation for this is that 
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VIMOS (A,, Map) 



VIMOS (Rv- Map) 



1" 



= 0.5 Ay 




= 2.5 Uy 



Figure 7. Maps of the CCM fits, the slope Ry (left) and normalization, Ay (right). The black rectangle is the aperture analyzed in H09 
(Figures 8 and 9). The red and blue circles denote the background and foreground galaxy nuclei, respectively to help with orientation. 
The VIMOS observations cover the dark spiral structures observed in the HST extinction map. The Ay values in both maps agree well. 
Realistic Ry values range from 0.8 to 5.2, on average much lower than the Milky Way value (Ry = 3.1). 



-25°07'34.0"- 




^0.5 Ay 



39.0"- 

21.80s 



21.60s 
RA 1J2000) 



0h48m21.40s 



Figure 8. The extinction (Ay) map based on the F606W HST 
image (black outline), using method B described in H09 (using a 
model background galaxy). Aperture II in the VIMOS analysis is 
marked with the green triangle. 



the spatial averaging is not over the whole HST aperture 
but sections of 0'.'67 squared. Unphysical fits (negative Ry) 
are found in those fibres closer to the foreground disk (fl4 
and fl5). The foreground disk contribution may not be non- 
negligible here (but we lack comparison fibres) and/or the 
spiral arm structure evident in the HST image may devi- 
ate from symmetry. The grayer extinction curve seen in the 



overlap region (Figure 7) is an observational effect of the 
VLT fibre's coarse spatial sampling. The seeing and IFU fi- 
bre scale (~ 0'.'7), each aperture has a linear size of ~ 0.9 
kpc in the foreground disk. At much higher spatial resolu- 
tion, we found a Milky Way-like extinction law using the 
Hubble data with ~ 60 pc resolution elements (e.g., H09). 
In contrast, the VLT fibre apertures each cover a much wider 
area, each of which include both optically thin and thick re- 
gions of the foreground disk (Figure 8). The disk opacity is 
dominated by the optical thick regions whose covering fac- 
tor does not change much with wavelength (see also Natta & 
Panagia 1984; Fischera & Dopita 2005, for the extinction law 
effects of an inhomogeneous screen). This resolution effect 
has been observed in other occulting pairs as well (Berlind 
et al. 1997; Keel & White 2001a,b). A Milky Way-type law 
occurs when the physical resolution is similar to the typical 
size of a Giant Molecular Cloud (^ 100 pc), still an order of 
magnitude below the VIMOS sampling. 

The Ay values (Figure 7, left) are in line with averages 
of the extinction maps presented in H09. Figure 8 shows 
the HST aperture with Aperture II overlaid. We note that 
the opaque structures in the HST extinction map (based on 
i^606M^ filter) dominate in Aperture II (e.g., the high value 
of Ay in fibres 6, 8, and 12). 

The fact that we still measure a slope in these fibres 
point to structures with similar extinction properties larger 
than what we found before (> 60 — 100 pc), but smaller 
than the VIMOS sampling (< 9 kpc), since Ry values are 
still lowered by oversampling. The spiral structures in the 
HST-hased extinction map (Figure 8) are a prime candidate 
for these structures. 

An alternative explanation is that the typical Ry values 
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Figure 9. Optical depth maps {Ay = 1.086 X ty) based on 
tlie three HST filters (F475W (B), F606W (V), and F814W (I). 
Optically thick structures stand out in the shorter wavelengths. 



of this galaxy's dusty ISM are higher than the canonical Rv 
to begin with (before resampling). In this case, the grain 
distribution responsible for the extinction curve would be 
distributed more towards smaller grains (more reddening), 
a sign of recently reprocessed or produced dust. If this were 
the case, this would point to a disk undergoing significant 



Measure 


Value 




Unit 


Note 


z 


0.064 








D 


263 




Mpc 




m - M 


37.1 








A^F475W 


-18.2 




mag. 


(a) 


A^F606W 


-18.9 




mag. 


(a) 


A^F814W 


-19.5 




mag. 


(a) 


B-V 


0.64 








V-I 


0.62 








log(M/Lv) 


0.26 




Lq/Mq 


(b) 


log{M/Lv) 


-0.6 




Lq/Mq 


(c) 


M, 


5.5 X 


109 


Mq 


(b) 


M, 


0.76 X 


109 


Mq 


(c) 


<M* > 


3.13 X 


109 


Mq 


(d) 


< r > 


0.122 






(d) 


n 


129.7 




arcsec^ 


(e) 


Mdust 


0.52 X 


106 


Mq 




Gas-to-dust 


~100-500 




(f) 


Mgas 


0.52- 


2.6 X 10** 


Mq 




/gas 


0.02- 


0.08 







Table 2. Characteristics of the Foreground Galaxy. 

(a) Measured from the HST image. 

(b) Bell & de Jong (2001), Figure 9, [B-V] - log{M/Lv)- 

(c) same but [V-I] - log(M/Lv) 

(d) Mean optical depth per pixel in the overlap region in H09. 

(e) Kron aperture angular size. 

(f) Gas-to-dust ratios vary across spiral disks (Smith et al. 2010; 
Leroy et al. 2011; Magrini et al. 2011) and these high values are 
typical for smaller disks. 

change, for example through the star-burst and shocks as- 
sociated with a close encounter. 



7 DISCUSSION 

With a new, more reliable, distance measurement for the 
foreground galaxy, we can use the HST observations to infer 
some of its global characteristics. Our aims were to gauge if 
the foreground galaxy is in any way unusual, specifically if 
it is interacting or if it is a fairly typical dwarf spiral which 
just happens to be positioned over a bright neighbour. 

We first estimate the stellar mass of the foreground 
galaxy. We calculate a V band absolute magnitude of My ~ 
— 18.9, based on a 9'.'0 (11.5 kpc) elliptical (Kron) aperture 
on the F606W HST irrmge and assuming the F606W is iden- 
tical to the Johnson V. Then, adopting mass-to-light ratios 
for V, based on either broadband B-V or V-I colour (Table 2) 
and the results of Bell &: de Jong (2001) in their Figure 9, a 
stellar mass of M* ~ 3.1 x IO^Mq. We can also estimate the 
likely circular velocity via the Tully-Fisher relation. Our V- 
band luminosity suggests Ujot ~ 97 km/ s (Verheijen 2001). 
All this suggest that the galaxy is a relatively low-mass spi- 
ral, comparable to M33 or NGC 4244 and consistent with 
the dust structure power spectrum observed. 

We now estimate the gas mass by first approximating 
the dust mass, and then assuming a dust-to-gas ratio. We 
derive the dust mass by assuming the mean optical depth we 
derived in H09 {jv ~ 0.12) is typical for the entire disk (i.e., 
the Kron aperture we used above). Implicitly, this assumes 
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that the dark features in front of the background galaxy 
are not tidal or asymmetric. The dust mass derived this 
way is a lower limit, as we assume a fiat distribution of a 
finite (truncated) disk, not an exponential one. The latter 
assumption is predicated on a scale-length of the dust larger 
than the stellar one, which appears to be valid for smaller 
disks (MacLachlan et al. 201f ; Holwerda et al. 2011; Holw- 
erda et al. 2012a). We then adopt a gas-to-dust ratio of 100 
and 500 (typical for lower mass disks in Leroy et al. 2011), 
yielding Mgas ~ 0.5 - 2.6 x 10* M©. 

Combining the stellar and gas masses, we infer a -rather 
uncertain- gas fraction /gas = Mgas/(Mgas -I- M,) ~ 0.02 — 
0.08. This gas fraction for a stellar mass of M, ~ IO^Mq 
is very low compared to typical values (Kannappan 2004; 
West et al. 2010), or the general trend of HI gas fraction 
with stellar mass (Catinella et al. 2010; Fabello et al. 2011, 
for bigger disks). However, because of the mass of dust avail- 
able for shielding within the Kron radius, the bulk of this 
gas is most likely in the molecular phase. Compared to the 
trend of molecular gas fraction with stellar mass (Saintonge 
et al. 2011), the gas fraction is still relatively low, although 
much less unusual. Assuming an atomic gas envelope of 
equal mass, in addition to the molecular disk within the 
Kron radius, the gas fraction would be typical for the mea- 
sured stellar mass. Spiral galaxies of similar stellar mass are 
in the HERACLES survey (Leroy et al. 2009), with similar 
molecular gas mass disks inferred from CO observations, all 
of which is well contained within the stellar radius (i?25). 

Recent Herschel results on the SED of stacked popula- 
tions of galaxies point to a low dust temperature (Tdust ~ 
15 K) and low optical/ultraviolet attenuation of the stel- 
lar population, together with a high relative dust mass 
(Afdust/Af*) from sub-mm emission for blue galaxies with 
a stellar mass of ~ 10® M© (Bourne et al. 2012). The dust 
geometry we observe in the foreground disk would explain 
these characteristics: dust is heated less by the stellar disk, 
because much of the dust is more radially extended that 
the stars. In addition, the stellar attenuation is not as high 
for the same reason, yet there is a relatively large amount 
of cold dust to emit at sub-mm wavelengths. This is remi- 
niscent of a nearby counterpart for this galaxy, NGC 4244 
(frot = 100 km /a), which also displays a large dust mass 
which is distributed in a clumpy medium throughout the 
disk, right up to the optical radius (i?25 MacLachlan et al. 
2011; Holwerda et al. 2011; Holwerda et al. 2012a). 

As discussed above, we suspect this pair is not inter- 
acting, although this scenario cannot be ruled out based on 
the respective redshifts. If the pair is physically closely as- 
sociated, the foreground galaxy's trajectory could take it 
through the halo of the more massive background galaxy. 
This, in turn, could have resulted in the extended dusty 
disk without the need for a strong tidal interaction. The in- 
teraction with the intergalactic medium and gaseous halo of 
the foreground galaxy would result in a displacement and re- 
processing of the ISM disk. However, if the pair is not inter- 
acting, the remaining interpretations are intriguing. What 
star formation history, in time and galactic radius, would 
result in such an extended dusty disk? Can radial migration 
(e.g., Roskar et al. 2010, Radburn-Smith et al. in prep) of 
the stars (inward) or ISM (outward), be the cause of the 
extend of the dusty disk? Or, very speculatively, is the stel- 
lar distribution in the core region physically thicker than 



the outer disk, and the whole galaxy disk is very dusty and 
opaque (a "heavy smoker" , see Disney et al. 1989; Phillips 
et al. 1991)? More detailed observations of the kinematics of 
this dusty disk (e.g., with ALMA or MeerKAT) would help 
settle these questions. 

The observed extinction curves can be well fit with a re- 
lation described in Cardelli et al. (1989) with grayer (lower 
Rv) slopes (Figure 7). We attribute this predominantly to 
the sampling of the VIMOS observations. In IFU observa- 
tions of a closer occulting pair that is interacting, UGC 
3995 A/B, we find a similarly gray extinction curve for sim- 
ilar physical sampling (~ 9 kpc, Holwerda et al. in prep). 
We note that in that pair, the distribution of Av values 
in the attenuation map based on the HST image display 
a much higher mean value than in this foreground spiral, 
pointing to a fundamental difference in the ISM disk. This 
difference in Av distribution between this pair and the in- 
teracting UGC 3995A/B, is circumstantial evidence that the 
foreground galaxy is not distorted by a gravitational inter- 
action with the background galaxy. Alternatively, the still 
high value of Rv measured over the course VIMOS sam- 
pling can be interpreted in tow ways: either more extended 
structures with the same attenuation properties exist in the 
disk (e.g., the spiral structures seen in Figure 8) or the native 
Rv values were much higher than the canonical Rv = 3.1, 
the signature of recently produced or reprocessed dust, cir- 
cumstantial evidence of recent interaction. 

Future work on occulting galaxy pairs will benefit 
greatly from higher spatial sampling of the next generation 
of IFU, e.g., the WEAVE instrument on the WHT or MUSE 
on VLT, both assisted with AO systems. To fully answer 
the question whether dust disks such as the one observed in 
the foreground disk is typical, a statistical sample of occult- 
ing galaxy pairs observed with an IFU instrument will be 
needed. 
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APPENDIX A: SPECTRAL SMOOTHING AND THE CCM FIT 

Table Al. The normalization {Ay) and slope {Rv) of the CCM fit to the attenuation curves in Aperture II, inferred from the raw and 
speetrally smoothed speetra (5,7,9,11 and 21 pixels respeetively). 



Fiber 


Raw 






13 A 






18.2 A 






23-4 A 






28-6 A 






54.6 A 






No. 


Av 


Rv 


X^ 


Av 


Rv 


X=^ 


Av 


Rv 


X^ 


Av 


Rv 


X = 


Av 


Rv 


X^ 


Av 


Rv 


X 




1 


0.59 


2.52 


2.07 


0.59 


2.53 


1.19 


0-59 


2-53 


0-91 


0-59 


2-53 


0-72 


0-59 


2.53 


0-59 


0.59 


2-53 


0-34 


2 


0.60 


1.71 


6.16 


0.60 


1.71 


3.36 


0.60 


1-71 


2-56 


0-60 


1-71 


2-05 


0-60 


1.71 


1-69 


0.60 


1-70 


0-85 


3 


0.67 


1.24 


24.48 


0.66 


1.24 


14.42 


0-66 


1-24 


11.17 


0-66 


1-23 


9-19 


0-66 


1-23 


7-87 


0.66 


1-21 


4-92 


4 


0.70 


2.96 


32.08 


0.69 


3.13 


14.92 


0-69 


3-18 


10.34 


0-69 


3-21 


7-83 


0-69 


3.23 


6.26 


0.69 


3-28 


3-00 


5 


0.52 


1.43 


37.72 


0.51 


1.48 


19.16 


0-51 


1-49 


13.51 


0-51 


1-50 


10-09 


0-51 
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Figure Al. The position of the fibres in Aperture II reported in Table Al and Figure B2. Grayscale was adjusted so clarify fibre labels. 
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Figure A2. Comparison between the raw and spectrally smoothed CCM fits, the slope (Rv) ^^^ normalization (Ay). 
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APPENDIX B: THE EFFECTS OF SPECTRAL SMOOTHING 
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Figure Bl. An example fibre from Aperture II, smoothed with different Hanning windows (13, 18.2, 23.4, 28.6, and 54.6 A) to illustrate 
the effect of smoothing the spectra to different resolutions. Our CCM fit from the raw fibres are robust in this Aperture and there is no 
need to smooth in the spectral direction to improve fit. 
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Figure B2. The raw attenuation curves for Aperture 11. The background (red) and the overlap (sohd black line) spectra. The observed 
attenuation curve (dashed line) and the fit to this line (dotted line). 
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